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Abstract 



We present an analysis of the neutrino oscillation solutions of the solar neu- 
trino problem in the framework of four-neutrino mixing where a sterile neu- 
trino is added to the three standard ones. We perform a fit to the full data 
set corresponding to the 825-day Super-Kamiokande data sample as well as 
to Chlorine, GALLEX and SAGE and Kamiokande experiments. In our anal- 
ysis we use all measured total event rates as well as all Super-Kamiokande 
data on the zenith angle dependence and the recoil electron energy spectrum. 
We consider both transitions via the Mikheyev-Smirnov-Wolfenstein (MSW) 
mechanism as well as oscillations in vacuum (just-so) and find the allowed 
solutions for different values of the additional mixing angles. This framework 
permits transitions into active or sterile neutrinos controlled by the addi- 
tional parameter cos^(t?23) cos^(t?24) and contains as limiting cases the pure 
fg-active and -sterile neutrino oscillations. We discuss the maximum al- 
lowed values of this additional mixing parameter for the different solutions. 
As a particularity, we also show that for MSW transitions there are solutions 
at 99 % CL at ■i?i2 mixing angles greater than j and that the best fit point 
for the zenith angle distribution is in the second octant. 
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I. INTRODUCTION 



Solar neutrinos were first detected already three decades ago in the Homestake experi- 
ment [Q and from the very beginning it was pointed out the puzzling issue of the deficit in 
the observed rate as compared to the theoretical expectation based on the standard solar 
model with the implicit assumption that neutrinos created in the solar interior reach 
the Earth unchanged, i.e. they are massless and have only standard properties and inter- 
actions. This discrepancy led to a change in the original goal of using solar neutrinos to 
probe the properties of the solar interior towards the study of the properties of the neutrino 
itself and it triggered an intense activity both theoretical as well as experimental, with new 
measurements being proposed in order to address the origin of the deficit. 

On the theoretical side, enormous progress has been done in the improvement of solar 
modeling and calculation of nuclear cross sections. For example, helioseismological obser- 
vations have now established that diffusion is occurring and by now most solar models 
incorporate the effects of helium and heavy element diffusion ^,0. iFiom the experimental 
point of view the situation is now much richer. Four additional experiments to the original 
Chlorine experiment at Homestake 1^ have also detected solar neutrinos: the radiochemical 
Gallium experiments on pp neutrinos, GALLEX and SAGE [0, and the water Cerenkov 
detectors Kamiokande and Super-Kamiokande [^|TU[. The latter have been able not 
only to confirm the original detection of solar neutrinos at lower rates than predicted by 
standard solar models, but also to demonstrate directly that the neutrinos come from the 
sun by showing that recoil electrons are scattered in the direction along the sun-earth axis. 
Moreover, they have also provided us with good information on the time dependence of the 
event rates during the day and night, as well as a measurement of the recoil electron energy 
spectrum. After 825 days of operation, Super-Kamiokande has also presented preliminary 
results on the seasonal variation of the neutrino event rates, an issue which will become 
important in discriminating the MSW scenario from the possibility of neutrino oscillations 
| Tl|JT^ . At the present stage, the quality of the experiments themselves and the 



m vacuum 



robustness of the theory give us confidence that in order to describe the data one must 
depart from the Standard Model (SM) of particle physics interactions by endowing neutri- 
nos with new properties. In theories beyond the SM, neutrinos may naturally have new 
properties, the most generic of which is the existence of mass. It is undeniable that the 
most popular explanation of the solar neutrino anomaly is in terms of neutrino masses and 



mixing leading to neutrino oscillations either in vacuum |13| or via the matter-enhanced 



MSW mechanism [|I^ . 

On the other hand, together with the results from the solar neutrino experiments we 
have two more evidences pointing out towards the existence of neutrino masses and mixing: 
the atmospheric neutrino data and the LSND results. The first one can be summarize in 
the existence of a long-standing anomaly between the predicted and observed U/^/ue ratio 
of the atmospheric neutrino fiuxes ||15[. In this respect it has been of crucial relevance 



the confirmation by the Super-Kamiokande collaboration |TO,ro of the atmospheric neu- 
trino zenith-angle-dependent deficit, which strongly indicates towards the existence of 
conversion. In addition to the solar and atmospheric neutrino results from underground 
experiments, there is also the indication for neutrino oscillations in the — >■ z/g channel by 



the LSND experiment [17|. All these experimental results can be accommodated in a single 
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neutrino oscillation framework only if there are at least three different scales of neutrino 
mass-squared differences. The simplest case of three independent mass-squared differences 
requires the existence of a light sterile neutrino, i.e. one whose interaction with standard 
model particles is much weaker than the SM weak interaction, so it does not affect the 



invisible Z decay width, precisely measured at LEP 1 18-22 



In this paper we present an analysis of the neutrino oscillation solutions of the solar 
neutrino problem in the framework of four-neutrino mixing where a sterile neutrino is added 
to the three standard ones. We perform a fit of the full data set corresponding to the 825- 
day Super-Kamiokande data sample as well as the data of the Chlorine, GALLEX and 
SAGE experiments. In our analysis we use all measured total event rates and all Super- 
Kamiokande data on the zenith angle dependence and the recoil electron energy spectrum. 
We consider both transitions via the Mikheyev-Smirnov-Wolfenstein (MSW) mechanism as 
well as oscillations in vacuum (just-so) and find the allowed solutions for different values of 
the additional mixing angles. Our analysis contains as limiting cases the pure z/e^active and 
z/g-sterile neutrino oscillations. We discuss the maximum allowed values of the additional 
mixing angles for which the different solutions are allowed. 

The outline of the paper is as follows. In Sec. |I| we summarize the main expressions for 
the neutrino oscillation formulas that we use in the analysis of solar neutrino data which take 
into account matter effects in the case of the MSW solution of the solar neutrino problem. 
We also present some improvement concerning the calculation of the regeneration of solar 
z/e's in the Earth in Appendix ^ Section |ITl| contains the summary of our calculations for 
the predictions of the different observables. Our quantitative results for the analysis of the 
four-neutrino oscillation parameters are given in Sec. |V[ Finally in Sec. |V| we summarize 
and discuss briefly our conclusions. 



II. FOUR NEUTRINO OSCILLATIONS 



In this paper we consider the two four-neutrino schemes that can accommodate the 
results of all neutrino oscillation experiments [^,^: 



Ami 



Ami 



(A) (nil < 1^2) > (ma < 7714) , (B) (mi < 1712) < (m^ < ml) 



(2.1) 



In both these mass spectra there are two pairs of close masses separated by a gap of about 
1 eV which gives the mass-squared difference AmggL = Am4^ responsible for the short- 
baseline (SBL) oscillations observed in the LSND experiment (we use the common notation 



Am2 . ; 



m^ — mj). We have ordered the masses in such a way that in both schemes 
Am^i produces solar neutrino oscillations and Am^tm = Am^g is responsible 



for atmospheric neutrino oscillations. With this convention, the data of solar neutrino 



experiments can be analyzed using the neutrino oscillation formalism presented in Ref. 



that takes into account matter effects. In this section we present the neutrino oscillation 
formulas that we use in the analysis of solar neutrino data. The transition probabilities 
that take into account matter effects in the case of the MSW solution of the solar neutrino 
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problem have been derived in Ref. [^]. Here we present some improvement concerning the 
calculation of the regeneration of solar u^s in the Earth (see the Appendix]^). 

In four-neutrino schemes the flavor neutrino fields u^l (« = e, s, /i, r) are related to the 
fields i>icL of neutrinos with masses m/c by the relation 

4 

UaL = '^Uak^kL (a = 6, S, jJ, t) , (2.2) 

k=l 

where [/ is a 4x4 unitary mixing matrix, for which we choose the parameterization 

U = UuU24U23UuUisUu, (2.3) 

where 

{Uij)ab = Sab + (cOST^ij - 1) {6iaSib + 6jaSjb) + Sm{}ij {6iaSjb - 6jakb) (2.4) 



represents a rotation in the i-j 2x2 sector by an angle ^ij. In the parameterization ( p.3| ) 
we have neglected, for simplicity, the possible presence of CP-violating phases. 

Since the negative results of the Bugey Pe disappearance experiment |Q imply that 
l^^esP + |t/e4p < 3 X lO'^ for Am^^L in the LSND-allowed region 0.2 eV^ < Am^g^ < 2 eV^ 
in the study of solar neutrino oscillations the matrices t/13 and f/i4 can be approximated 
with the unit matrix [i.e. di^ = = 0) and we obtain 



U = U^^U2,U2zU^2. (2.5) 



Explicitly, we have 



U 



( Ci2 S12 

-S12C23C24 C12C23C24 S23C24 S24 

S12(C23«24S34+S23C34) -Cl2(s23C34+C23«24S34) C23C34-S23.524S34 C24S34 

\ S12(C23«24C34-S23^*34) C12 (s23.534-C23.524C34) -(C23.534+S23.524C34) C24C34 } 



(2.6) 



where -Qyi, ^23, ^24^, ^^34 are four mixing angles and c^- = cos ^ij and Sij = sin 

Since solar neutrino oscillations are generated by the mass-square difference between 
z/2 and z/i, it is clear from Eq. ( |2.6| ) that the survival of solar z/g's mainly depends on the 
mixing angle ■&12, whereas the mixing angles ■}}23 and -(924 determine the relative amount of 
transitions into sterile or active and u^-. Let us remind the reader that and Pj- 
cannot be distinguished in solar neutrino experiments, because their matter potential and 
their interaction in the detectors are equal, due only to neutral- current weak interactions. 
The active/sterile ratio and solar neutrino oscillations in general do not depend on the mixing 
angle 'd^^, that contribute only to the different mixings of and z/^, and depends on the 
mixing angles -^23 "^24 only through the combination cos'i923 cos'i924- Indeed, from Eq. ( p.6|) 
one can see that the mixing of Ug with ui and z/2 depends only on 1912 and the product 
cos ^^23 cost? 24- Moreover, instead of and Ur, one can consider the linear combinations 




— sin -i? — cos f sin -(934 cos 'd. 
cos — sin / I cos 1^)34 — sin -^s^ 




(2.7) 



4 



with 



tan^ = 2.8 



The mixing of Ua and Ub with z/i and 1^2 is given by 



Ual = -S12 Y 1 - , Ua2 =012^1- 4^4^^ , Ubl = f/fe2 = . (2.9) 

Therefore, the oscillations of solar neutrinos depend only on i^u and the product 
cos'i923 cos'i924- If COS '(923 COS "(924 7^ 1) solar UgS can transform in the linear combination 
Ua of active z/^ and u^-. We distinguish the following limiting cases: 

• If cos'(923 cos -(9 24 = then Usi = Us2 = 0, Uai = — sin'i9i2, Ua2 = cos'du, corresponding 
to the limit of pure two-generation z/g z/^ transitions. 

• If cos "(923 cos ^^24 = 1 then Usi = — sin'i9i2, Us2 = cos'i9i2 and Uai = Ua2 = and we 
have the limit of pure two-generation Uf, transitions. 

Since the mixing of Ve with vi and z/2 is equal to the one in the case of two- generations 
(with the mixing angle 'di2)-, the mixing of Vg with vi and z/2 is equal to the one in the case 
of two-generations times cos "(923 cos d2i and the mixing of z/^ with z/i and V2 is equal to the 
one in the case of two-generations times \/1 — cos^ '§2^ cos^ ^924, it is clear that in the general 
case of simultaneous and z/g — > Ua oscillations the corresponding probabilities are 

given by 

P^:Z.^ = cos' ^23 COS^ .924 (l - P^^Xj , (2-10) 
P^:X. = (1-C0S2^23C0S2^24) {l ' P^^Zj • (2.11) 

These expressions satisfy the relation of probability conservation P^^^^^ + Pu^Zus + 

pS\in _ 



Ve^Va 



If Am^i is in the MSW region (10 ^eV^ < Am2i < 3 x 10 ^eV^), the survival proba- 
bilities of solar z/e's is given by p2[ 



p!:Zu. = ^ + Q - p. ) cos 2^12 cos 2^f^ . (2.12) 

Here the angle is the effective mixing angle in matter corresponding to the vacuum 
mixing angle di2 and given by 



tan 2^f^ = , / , (2. 13) 

1 - /l/AmiiCos279i2 

with 

A = Acc + cos' cos' Anc ■ (2.14) 
The quantities Acc and Aj^c describe the matter effects and are given by 
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Acc = 2V2GfEN, , Anc = -V2GFENn , 



(2.15) 



where and iV„ are, respectively, the number densities of electrons and neutrons in the 
medium, E is the neutrino energy and Gp is the Fermi constant. The effective mixing angle 
'd^2 Eqs. (|2.12| ) and ( p.lOj ) must be evaluated at the point of neutrino production inside 
of the Sun. The quantity P^. in Eq. ( |2.12| ) is the crossing probability given by the usual 
two-generation formula (see [^) and the replacement of the two-generation expression for 
A with that given in Eq. ( ^.141 ). 

During the night solar neutrinos cross the Earth before reaching the detector and re- 
generation of z/e's is possible 



23 . In the four-neutrino schemes under consideration, the 



probabilities of 



and z/p 



transitions after crossing the Earth are given by p2| 



P: 



Sun+Earth 



Sun 



+ 



(1 - 2P.^:i:.j (p, 



Earth 



sm' 



COs2't?i2 



(2.16) 



pSun+Earth ^ pSun , (2^.^-.. " COS^ ^,3 COS^ ^,4) (P.^-^t ' ^O^' ^12 COS^ ^23 COS^ .^,4) 
^,^us^ COs2i9i2COs2i923COs2^?24 

(2.17) 

The probability of ~^ ^'a transitions is given by the conservation of probability: 

pSun+Earth _ 1 pSun+Earth pSun+Earth 

The probabilities P^^^^^ and P^^^^^ in Eqs. ( p.l6|) and ( |2.17| ) can be calculated by inte- 



grating numerically the differential equation that describes the evolution of neutrino flavors 
in the Earth (see p2[) or by using the analytical solution assuming a step-function profile of 



the Earth matter density (see ||2^) . However, we notice that the probabilities P^^^I^^^ and 
Pu2^^us independent, because, as shown in the Appendix 0, they are related by 

= COS^ ^23 COS^ ^24 (1 - PS'^t) . (2.18) 

Therefore, in the analysis of solar neutrino data we need to calculate only P^^^_^^^. 

If Am2i is in the range of the vacuum oscillation solution of the solar neutrino problem 
(lO^^^eV^ < Amli < 10^^ eV^), the survival probability of solar is given by the two- 
generation formula 

^:"-:.. = l-sin^2^,2sin^^^, (2.19) 

where E is the neutrino energy and L is the Sun-Earth distance, whose seasonal variations 
must be taken into account. In this case there is no matter effect during neutrino propagation 
in the Earth. 



III. DATA AND TECHNIQUES 

In order to study the possible values of neutrino masses and mixing for the oscillation 
solution of the solar neutrino problem, we have used data on the total event rates measured 
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in the Chlorine experiment at Homestake 0, in the two Galhum experiments GALLEX 
and SAGE [^,0 and in the water Cerenkov detectors Kamiokande and Super-Kamiokande 
shown in Table |. Apart from the total event rates, we have in this last case the zenith 
angle distribution of the events and the electron recoil energy spectrum, all measured with 
their recent 825-day data sample W^. Although, as discuss in Ref. p6[ the inclusion of 



|10|] . Although, as discuss in Ref. 
Kamiokande results does not affect the shape of the regions, because of the much larger 
precision of the Super-Kamiokande measurement, it is convenient to introduce it as in this 
way the number of degrees of freedom for the fit of the rates only is 4 — 3 = 1 (instead of zero 
degrees of freedom), that allows the construction of a well-defined Xmin confidence level. 

For the calculation of the theoretical expectations we use the BP98 standard solar model 
of Ref. 1^ . The general expression of the expected event rate in the presence of oscillations 



in experiment i in the four-neutrino framework is given by K. 



th 



k=l,8 

+cr.,(E,)(l-(P,,^,,)-(P,,^,j) 



(3.1) 



where Ei, is the neutrino energy, (pk is the total neutrino flux and is the neutrino energy 
spectrum (normalized to 1) from the solar nuclear reaction k with the normalization given 
in Ref. Here a^^i (o"a;,i) is the z/g {ux, x = fi, t) interaction cross section in the Standard 
Model with the target corresponding to experiment i. For the Chlorine and Gallium experi- 
ments we use improved cross sections aa,i{E) {a = e, x) from Ref. [^. For the Kamiokande 



and Super-Kamiokande experiment we calculate the expected signal with the corrected cross 
section as explained below. {Piy^_,^^) is the time-averaged survival probability. In case of 
MSW transitions Pu^^ue Pv^^us ^"^^ given in Eqs.( p.l(^ ) and (|2.17| ) respectively. 

For vacuum oscillations we must include the effect of the Earth orbit eccentricity. The 
yearly averaged probabihty is obtained by averaging Eg. (|2.19|) with L{t) = Lo[l — e cos 27i-^]: 



(P.,^.J = (P,s- ) = 2-sin^2^ 



12 



1 



COS 



2E 



Jo 



eAml^LQ 
2E 



(3.2) 



where s is the orbit eccentricity (0.0167), Lq is the average Earth orbit radius (1.496 x 10^ 
km) and Jo{x) is the Bessel function. We have also included in the fit the experimental 
results from the Super-Kamiokande Collaboration on the zenith angle distribution of events 
taken on 5 night periods and the day averaged value, which we graphically reduced from 
Ref. [|lOl- For MSW oscillations we compute the expected event rate in the period a in the 



presence of oscillations as, 

-j^ l'T{cOS^rnax.a) 

Rlia = / drY.'Pk jdE, Xk{E,) X (J,^sk{E,){P,,^,M) 



^'^a JricOS^rnin.a) k = l 8 

+axME,)(l - (P..^..(r)) - (P..^..(r)) 



(3.3) 



where r measures the yearly averaged length of the period a normalized to 1, so Ar^ = 
T{cos^rnax,a) — t(cos $mm,a) = .500, .086, .091, .113, .111, .099 for the day and five night 
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periods. Notice that for vacuum oscillations there is no matter effect during neutrino propa- 
gation in the Earth. In this case Rl'l^ — -^sfc> given in Eq. (|3.1| ). The Super-Kamiokande 
Collaboration has also presented the results on the day-night variation in the form of a day- 
night asymmetry. Since the information included in the zenith angle dependence already 
contains the day-night asymmetry, we have not added the asjTumetry as an independent 
observable in our fit. 

The Super-Kamiokande Collaboration has also measured the recoil electron energy spec- 
trum. In their published analysis after 504 days of operation they present their results 
for energies above 6.5 MeV using the Low Energy (LE) analysis in which the recoil en- 
ergy spectrum is divided into 16 bins, 15 bins of 0.5 MeV energy width and the last bin 
containing all events with energy in the range 14 MeV to 20 MeV. Below 6.5 MeV the back- 
ground of the LE analysis increases very fast as the energy decreases. Super-Kamiokande 
has designed a new Super Low Energy (SLE) analysis in order to reject this background 
more efficiently so as to be able to lower their threshold down to 5.5 MeV. In their 825-day 
data |jlO| they have used the SLE method and they present results for two additional bins 
with energies between 5.5 MeV and 6.5 MeV. In our study we use the experimental results 
from the Super-Kamiokande Collaboration on the recoil electron spectrum divided in 18 
energy bins, including the results from the LE analysis for the 16 bins above 6.5 MeV and 
the results from the SLE analysis for the two low energy bins below 6.5 MeV. The general 
expression of the expected rate in a bin in the presence of oscillations, R^^, is similar to that 
in Eq. (|3.1|) , with the substitution of the cross sections with the corresponding differential 
cross sections folded with the finite energy resolution function of the detector and integrated 
over the electron recoil energy interval of the bin, T^in <T< T^a.x- 



The resolution function Res{T, T') is of the form |^,^ 

1 



dV 



Res{T, T') 



27r(0.47VT'(MeV)) 



exp 



(T-T 



l\2 



0.44T'(MeV)_ 



(3.4) 



(3.5) 



and we take the differential cross section daa{Ejj, T')/dT' from pO[ . 

In the statistical treatment of all these data we perform a analysis for the different 
sets of data, following closely the analysis of Ref. with the updated uncertainties given 
in Refs. |]26|-p8|, as discussed in Ref. [^]. We thus define a function for the three set of 
observables xLes, xLnith, and Xspectrum where in both xLith, and Xspectrum we allow for a 
free normalization in order to avoid double-counting with the data on the total event rate 
which is already included in Xratos- combinations of observables we define the of 

the combination as the sum of the different x^'s- In principle such analysis should be taken 
with a grain of salt as these pieces of information are not fully independent; in fact, they 
are just different projections of the double differential spectrum of events as a function of 
time and energy. Thus, in our combination we are neglecting possible correlations between 
the uncertainties in the energy and time dependence of the event rates. 
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IV. RESULTS 



As explained in Sec. |I|, for the mass scales invoked in the explanation of the atmospheric 



and LSND data and after imposing the strong constraints from the Bugey [33| and CHOOZ 



3J] reactor experiments, the relevant parameter space for solar neutrino oscillations in the 
framework of four-neutrino mixing is a three dimensional space in the variables Am^i, ^9i2 
and cos^('i923) cos^('(924) = cl^c^^- As shown in Section |T|, the case cl^c^^ = corresponds to 
the usual two-neutrino oscillations z/g z/^ where Ua is the admixture of and given 
in Eq. ( p.7|) , thus an active neutrino. The other extreme case cl^c^^ = 1 corresponds to the 
usual two-neutrino oscillations of Uf, into a pure sterile neutrino. 

In our choice of ordering the neutrino masses in the two schemes (|2.1|) the mass-squared 
difference Am^i is positive. The mixing angle -^u can vary in the interval < '&12 < |. In 
the case of vacuum oscillations, the transition probabilities are symmetric under the change 
^12 ^ f ~ ''^12 and each allowed value of sin^(2'(9i2) corresponds to two allowed values of 
■&12. On the other hand, in the case of the MSW solutions the transition probabilities are 
not invariant under the change i!}i2 ^ f — ^12 and resonant transitions are possible only for 
values of 'di2 smaller than |. In the analysis of the observable rates, we present the results 
in the common plot of sin^(2'i9i2) due to the fact that the allowed region does not extend 
to '(?i2 > J- But, when we include the rest of observables we remark that this is not the 
case, and we present the results as a function of sin^('(9i2) (see ||3^) showing that there is 



a portion of the space of parameters in the second octant of •^}i2 allowed at 99 % CL. This 
enlarged parameter space is also used in Ref. p^ . 

We first present the results of the allowed regions in the three-parameter space for the 
different combination of observables. In building these regions, for a given set of observables, 
we compute for any point in the parameter space of four-neutrino oscillations the expected 
values of the observables and with those and the corresponding uncertainties we construct 
the function x^(Am^25 ^12, C23C24)obs- We find its minimum in the full three-dimensional 
space considering as a unique framework both MSW and vacuum oscillations. The allowed 
regions for a given CL are then defined as the set of points satisfying the condition: 

x\Ami„ ^12, 44),,, - xLn,ots < Ax'(CL, 3 dof) (4.1) 

where, for instance, Ax^(CL, 3 dof)=6.25, 7.83, and 11.36 for CL=90, 95, and 99 % respec- 
tively. In Figs. |l|-0 we plot the sections of such volume in the plane (Am2]^, sin^(2't9i2)) or 
(Am^i, sin^('t9i2)) for different values of 44. 

Figures ^ and ^ show the results of the fit to the observed total rates only. We find that 
both at 90 and 99% CL, the three-dimensional allowed volume is composed of three sepa- 
rated three-dimensional regions in the MSW sector of the parameter space (Fig. P, which 
we denote as SMA, LMA and LOW solutions following the usual two-neutrino oscillation 
picture and a "tower" of regions in the vacuum oscillations sector (Fig. H). The values of 
the minimum of the in the different regions are given in Table |T[ The global minimum 
used in the construction of the volumes lies in the SMA region and for a non-vanishing 
value of C23C24 = 0.3, although, as can be seen in the first panel in Fig. ||, this is of very 
little statistical significance as A^^ for the SMA solution is very mildly dependent on C23C24 
(Ax2<0.5for 44 < 0.5). 



9 



As seen in Fig. |l|, the SMA region is always a valid solution for any value of C23CI4. 
This is expected as in the two-neutrino oscillation picture this solution holds both for pure 
active-active and pure active-sterile oscillations^. On the other hand, both the LMA and 
LOW solutions disappear for a value of the mixing c^^cl^ > 0.5(0.3). Unlike active neutrinos 
which lead to events in the water Cerenkov detectors by interacting via neutral current with 
the electrons, sterile neutrinos do not contribute to the Kamiokande and Super-Kamiokande 
event rates. Therefore a larger survival probability for neutrinos is needed to accommo- 
date the measured rate. As a consequence a larger contribution from neutrinos to the 
Chlorine and Gallium experiments is expected, so that the small measured rate in Chlorine 
can only be accommodated if no ^Be neutrinos are present in the flux. This is only possible 
in the SMA solution region, since in the LMA and LOW regions the suppression of ^Be 
neutrinos is not enough. 

In Table fT| we give the maximum values of C23C24 for which the different solutions are 
allowed at the 90 and 99% CL according to different statistical criteria which we discuss 
below. In Fig. ^ we plot the corresponding sections in the vacuum oscillation sector. As 
seen in the figure, as C23C24 grows, the vacuum oscillation solution becomes more restricted 
in the allowed values of mass splittings till becoming a narrow band at Amg^ ~ 10~^° eV^ 
for the pure sterile case. 

Figure ^ shows the regions allowed by the fit of both total rates and the Super- 
Kamiokande zenith angular distribution in the MSW sector of the parameter space. In 
the vacuum oscillation section no day-night variation is expected. Also plotted is the ex- 
cluded region at 99 % CL from the zenith angular measurement. This exclusion volume is 
built as the corresponding three-degree-of-freedom region for the of the zenith angular 
data with respect to the minimum value Xmin zen = 0-8 which occurs at Am^i = 2.7 x 10~^ 
eV^, sin^('(9i2) = 0.85, and C23C24 = 0.0. We remark that this minimum is placed in the 
second octant and this was not included in past analysis of two-fiavor MSW solutions al- 
though it leads to little effect in the final results of the allowed regions. As seen in the figure 
and also in Table |nT| , the main effect of the inclusion of the day-night variation data is to 
cut down the lower part of the LMA region and to push towards slightly higher values the 
maximum CI3C24 for which the LMA and the LOW solutions are still valid. 

In Figs. § and ^ we plot the regions allowed by the fit of both total rates and the Super- 
Kamiokande energy spectrum. Also plotted is the excluded region at 99 % CL from the 
spectrum data which is obtained as the corresponding three-degree-of-freedom region for 
the of ths spectrum data with respect to the minimum value Xmin spec = 15.1 which 
occurs in the vacuum solution sector at Am^i = 6.3 x 10~^° eV^ and sin^(2'(9i2) = 1, and it 
is almost independent of C23CI4. As seen in the figure and also in Table [m| , the main effect 



Notice, however, that the statistical analysis is different: in the two-neutrino picture the pure 
active-active and active-sterile cases are analyzed separately, whereas in the four-neutrino picture 
they are taken into account simultaneously in a consistent scheme that allows to calculate the 
allowed regions with the prescription given in Eq. ( |4.1| ). We think that the agreement between the 
results of the analyses with two and four neutrinos indicate that the physical conclusions are quite 
robust. 
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of the inclusion of the spectrum data in the MSW regions is also to push towards shghtly 
higher values the maximum C23C24 for which the LMA and the LOW solutions are still valid. 
Figure ^ shows that the LMA region at 99 % CL extends to high values of Am^i , even above 
10~^eV^ for C23CI4 < 0.1. Since the atmospheric mass squared difference Aml^^ lies between 
10^'^ and 10~^eV^ (see [0), one may wonder if the solar and atmospheric mass squared 
differences may coincide and three massive neutrinos may be enough for the explanation 
of solar, atmospheric and LSND data. The answer to this question is negative, because 
in the high-Am^i part of the 99 % CL LMA region the mixing angle 621 is large, 0.3 < 
sin^(6'2i) < 0.7, and in this case disappearance of z/g's should be observed in long-baseline 
reactor experiments, contrary to results of the CHOOZ |Q experiment. In other words, 
the results of the CHOOZ experiment, that have not been taken into account in the present 
analysis, forbid the part of the 99 % CL LMA region that extends above Am^i — 10"^ eV^. 
For this reason we cut the plots at this value. For the vacuum sector, once the spectrum 
data is included the higher Arn^i are favored but we find no region at the 90% CL for any 
value of CI3C24 > 0.2 and at the 99% CL the region totally disappears for C23C24 > 0.8. 

Figures ^ and |^ show the results from the global fit of the full set of data. The values of 
the minimum iii the different regions are given in Table The global minimum used in 
the construction of the volumes lies in the LMA region and for vanishing c\^c^^ corresponding 
to pure z/g-active neutrino oscillations. 



In Table |T| we give the maximum values of c\^c\^ for which the different solutions are 



allowed at the 90 and 99% CL according to different statistical criteria. The use of each 
criteria depends on the physics scenario to which the result of our analysis is to be applied. 

• Criterion 1 (CI): The maximum allowed value of the mixing c\^c\^ at a given CL 
for a given solution is defined as the value for which the corresponding region of the 
allowed three-dimensional volume defined as a 3-d. o.f shift with respect to the global 
minimum in the full parameter space, disappears. In the first row in Fig. ^ we plot 
the values of Ax^ defined in this way for the different solutions as a function of C23C24 
for the fit of the total rates only (left panels) and for the global analysis (right panels). 
This criterion is the one used in building the regions in Figs. |1]-|^. It is applicable to 
models where no region of the parameter space MSW-SMA, MSW-LMA, MSW-LOW 
or vacuum is favored. 

• Criterion 2 (C2): The maximum allowed value of the mixing c\^c\^ at a given CL 
for a given solution is defined as the value for which the corresponding allowed three 
dimensional region defined as a 3-d. o.f shift with respect to the local minimum in the 
corresponding region, disappears. In the second row in Fig. |^ we plot the values of Ax^ 
defined in this way for the different solutions as a function of C23C24 for the fit to the 
total rates only (left panels) and for the global analysis (right panels). This criterion 
holds for models where only a certain solution, MSW-SMA, MSW-LMA, MSW-LOW 
or vacuum is possible and it yields less restrictive limits. 

• Criterion 3 (C3): The maximum allowed value of the mixing c\^c\^ at a given CL is 
obtained calculating the two-dimensional allowed regions in the sin^ ^i2-Am^2 plane 
for each fixed value of C23CI4. These allowed regions are defined through the 2-d.o.f shift 
with respect to the global minimum in the plane sin^ ^i2-Am^2 (this is the analogous 
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of Criterion 1 for two parameters). For each solution, the maximum allowed value of 
C23C24 is that for which the corresponding two dimensional region in the sin^ ^i2-Amf2 
disappears. This criterion is the equivalent to the usual two-neutrino analysis but with 
z/g oscillating into a state which is a given superposition of active and sterile neutrino. 

In the third row in Fig. ^ we plot the difference Ax^ between the local minimum of 
for each solution and the global minimum in the plane sin^ ^i2-Am^2 ^ function 

of 

^23^24- Notice that for the analysis of rates only the minimum in the plane sin^ 612 
and Am^2 occurs always in the MSW-SMA region for any value of c\^c\^- Thefore the 
curve for the MSW-SMA solution corresponds to the horizontal Ax^ = line and it is 
not shown. For the global analysis, when C23CI4 < 0.1 the minimum in the plane occurs 
for the MSW-LMA (dashed line) solution while for C23C24 > 0.1 it moves to the MSW- 
SMA (full line). For this reason the curve for the MSW-SMA (MSW-LMA) solution 
is only seen for C23CI4 < 0.1 (> 0.1) while for c\^C2^ > 0.1 (< 0.1) it coincides with the 
Ax^ = line. In general this criterion is applicable for models where the additional 
mixing (?2->(?2^ fixed a priori to some value, so that the model in fact contains only 
two free parameters (if C23C24 is larger than the maximum allowed for a given solution, 
it means that that solution is not allowed in the specific model). 

V. SUMMARY AND DISCUSSION 

At present, the Standard Model assumption of massless neutrinos is under question due 
to the important results of underground experiments. Altogether they provide solid evidence 
for the existence of anomalies in the solar and atmospheric neutrino fiuxes which could be 
accounted for in terms of neutrino oscillations z/g — and Vp,^ Vx-, respectivelyly. Together 
with these results there is also the indication of neutrino oscillations in the Vf. channel 

obtained in the LSND experiment. All these experimental results can be accommodated 
in a single neutrino oscillation framework only if there are at least three different scales of 
neutrino mass-squared differences. The simplest way to open the possibility of incorporating 
the LSND scale to the solar and atmospheric neutrino scales is to invoke a sterile neutrino, 
i.e. one whose interaction with Standard Model particles is much weaker than the SM weak 
interaction, so that it does not affect the invisible Z decay width, precisely measured at 
LEP. The sterile neutrino must also be light enough in order to participate in the oscilla- 
tions involving the three active neutrinos. After imposing the present constrains from the 
negative searches at accelerator and reactor neutrino oscillation experiments one is left with 
two possible mass patterns which can be included in a single four-neutrino framework as 
described in Sec. 0. 

In this paper we have performed an analysis of the neutrino oscillation solutions to 
the solar neutrino problem in the framework of four-neutrino mixing. We consider both 
transitions via the Mikheyev-Smirnov-Wolfenstein (MSW) mechanism as well as oscillations 
in vacuum. In what solar neutrinos is concerned our formalism contains one additional 
parameter as compared to the pure two-neutrino case: cos^('(?23) cos^('(924), where 'Q2J, and 
•Q2A give the projections of the sterile neutrino into each of the two heavier states responsible 
for explanation to the atmospheric neutrino anomaly. In this way, the formalism permits 
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transitions into active or sterile neutrino controlled by the additional parameter and contains 
as limiting cases the pure z/g-active and z/g -sterile neutrino oscillations. 

We have studied the evolution of the different solutions to the solar neutrino problem in 
this three parameter space when the different set of observables are included. In Figs. 
we plot the sections of such volume in the plane (Am2i, sin^(2'(9i2)) or (Am2i, sin^('i9i2)) for 
different values of C23C24. As a particularity, we also show that for MSW transitions there 
are solutions at 99 % CL at mixing angles greater than j and that the best fit point for 
the zenith angle distribution is in the second octant. 

Our results show that the SMA region is always a valid solution for any value of C23C24. 
This is expected as in the two-neutrino oscillation picture this solution holds both for pure 
active-active and pure active-sterile oscillations. On the other hand, the LMA, LOW and 
vacuum solutions become worse as the additional mixing C23C24 grows and they get to dis- 
appear for large values of the mixing. The main quantitative results of our analysis are 
summarized in Table fT| and Fig. |^ where we give the maximum values of C23C24 for which 
the different solutions are allowed at the 90 and 99% CL according to different statistical 
criteria which depend on the physics scenario to which the result of our analysis is to be 
applied. 
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APPENDIX A: DERIVATION OF THE RELATION BETWEEN P^^Z^}" AND 

pEarth 

In this appendix we derive the relation ( |2.18| ) between P^^^^^ and P^^^_^^^. In order to 
see the reason of this relation, it is useful to consider the most general 4x4 mixing matrix 
(without CP-violating phases) given in Eq. ( p.3|) , that can be written as 

U = U'U,2 , (Al) 

with 

U' = UuU24U23UuU,3. (A2) 

Let us define the neutrino states 

\K)= E ^-1^") (^ = 1,2,3,4). (A3) 

The amplitudes of transitions in the Earth for /c = 1, 2 are given by 

4 4 

= {^ASW) = Y.^i^M){K\sW) = J2uLrS'r,, (A4) 

r=l r=l 

where the unitary operator S describes the evolution inside the Earth and S'^.). = (i^^|5|z/fc). 

It has been shown in Ref. |^ that the matter effects inside the Earth can generate only 
transitions between ui, V2 and i^J, z^2- Then, we have 

^^T-^Vcc ^ ^al ^Ik + U'a2 S'2k = 1; 2) , (A5) 

= = U^k (A; = 3, 4), (A6) 



and the transition probabilities are given by 

P^^^-^Ua = ^al I'S'lfcP + l^'Zkl'^ + 2f^Qi '^^[Slk^2k\ = 2) , (A7) 

P^:^t=U'^k ik = 3A). (A8) 

Since the evolution operator S is unitary, the matrix 5" is unitary and we have the 
relations 

I C' 1 2 _ I c' 1 2 — p' I C' 1 2 _ I c' 1 2 _ 1 p/ Ql C'* I Q' Q'* _ fi /' A Q^ 

1^121 ~ l^2ll —-'12 5 I*-"!!! ~ 1^221 ~^ -'12 J *J21 "r ^12 ^22 ~ ' l^^J 

where is the probability of 1^2 ^ ly'i transitions, that is equal to the probability of 
^ V2 transitions. It is easy to check that the unitarity constraints (|A9|) are equivalent to 
the probability conservation relations 

EC".l = i' p^^ = i. (Aio) 

fc=l a=e,s,fi,T 
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Let us notice that by construction the matrix U' is such that 

[7^2 = 0, (All) 

and the probabihty of Vk transitions inside the Earth depend only on f/^^ = 

cos t^f 3 cos -(9^4 and P'y^- 

P^,Zl = U'A (1 - P[,) P^-^l = U'l P[, . (A12) 

On the other hand, in general, for a given mixing matrix f/, the transition probabilities 
Pf^^^^ with a = s,fi,T depend on two independent quantities, A^g['S'n'S'2i]. Hence, 

in general the probabihties -P^^*,!?^ and P^,^^^^ must be calculated independently. However, 
if 

KiK2 = {a = s,fi,T), (A13) 
also P^^f^^^ depends only on the elements of the mixing matrix and P^j- 

P^:^l = K\ (1 - + Pu , C"'. = UL\ Pu + Kl (1 - Pu) , (AM) 

for a = s, fi, T. 

In the approximation i^ia = i)i4 = that we use in the analysis of solar neutrino data, 
we have 

K, = l K, = {a = s,^I,T). (A15) 

Therefore, the condition ( [A13| ) is satisfied and we obtain 

pEarth _ i p/ pEarth _ p' ^'A1R'\ 



P^:^i = U'^2 Pu , = U'^2 (1 - Pl2) ■ (A17) 

Eliminating P^j from the relations ( |A16| ) and ( [AIT] ), we obtain 

P^:-^t = Kl (1 - {k = l,2;a = s, /i, r) . (A18) 

In particular, for k = 2 and a = s, we obtain the useful relation ( |2.18| ) between P^^*^ 
and P^;Z%- 

In general, considering the possibility of small but non-zero i^ia and/or i^i^, if the mixing 
angles are such that U'^^ = for a = s, /i, r, we have the relation 

(pEarth \ 
1-^1^ J (« = ^,/x,r), (A19) 

whereas if = we obtain the relation 



PuX. = Sp":-^". (« = r) . (A20) 



jEarth ^ al pEarth 
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Experiment 


Rate 


Ref. 


Units 


^BP98 


Homestake 
GALLEX + SAGE 
Kamiokande 
Super-Kamiokande 


2.56 ±0.23 
72.3 ±5.6 
2.80 ±0.38 
2.45 ± 0.08 


i 

m 

i 

[10| 


SNU 
SNU 

10^ cm-2 
10^ cm-2 g-i 


7.8 ± 1.1 
130 ±7 
5.2 ±0.9 
5.2 ±0.9 



TABLE I. Measured rates for the Chlorine, Gahium, Kamiokande and Super-Kamiokande 
experiments. 







rates 


rates±zenith 


rates±spectrum 


global 






5.6 X 10-^ 


6.2 X 10-^ 


5.2 X 10-6 


5.2 X 10-6 


SMA 


sin2(2^9) 


6.0 X 10-3 


5.4 X 10-3 


4.7 X 10-3 


4.7 X 10-3 




„2 2 
•^23^24 


0.3 


0.4 


0.0 


0.0 




xLin (Prob %) 


0.55 (45) 


6.2 (40) 


23.9 (16) 


29.7 (16) 






1.4 X 10-^ 


4.2 X 10-^ 


1.4 X 10-^ 


3.9 X 10-^ 


LMA 


sin2(2i?) 


0.68 


0.81 


0.68 


0.78 




„2 2 
'-23^24 


0.0 


0.0 


0.0 


0.0 




xLn (Prob %) 


3.80 (5) 


8.6 (20) 


23.1 (19) 


29.1 (18) 




Am^ 


1.3 X lO-'^ 


1.1 X lO-'^ 


1.0 X 10-^ 


1.0 X 10-^ 


LOW 


sin2(2i?) 


0.93 


0.93 


0.93 


0.93 




„2 2 
''23''24 


0.0 


0.0 


0.0 


0.0 




xLn (Prob %) 


8.3 (0.4) 


13.6 (3.4) 


27.9 (6.4) 


33.0 (8.1) 




Am^ 


9.1 X 10-^° 


9.1 X 10-^° 


4.5 X 10-^° 


4.4 X 10-^° 


Vacuum 


sin2(2T?) 


0.78 


0.78 


0.9 


0.9 




„2 2 
•-23^24 


0.0 


0.0 


0.0 


0.0 




X^„, (Prob %) 


4.5 (3.4) 


9.9 (13) 


28.8 (5.1) 


34.3 (6.1) 



TABLE IL Best fit points and the corresponding probabihties for the different solutions to the 



solar neutrino deficit and for different combinations of observables. 
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Observ. 


Grit. 


LMA 90 (99) 


LOW 90 (99) 


VAG 90 (99) 




Gl 


0.25 (0.54) 


- (0.36) 


0.56 (1.0) 


Rates 


G2 


0.44 (0.69) 


0.54 (1.00) 


1.0 (1.0) 




G3 


0.14 (0.43) 


- (0.18) 


0.12 (1.0) 




Gl 


0.30 (0.61) 


- (0.45) 


0.61 (1.0) 


Rates+SK Zenith 


G2 


0.45 (0.74) 


0.59 (0.90) 


1.0 (1.0) 




G3 


0.19 (0.49) 


- (0.23) 


0.15 (1.0) 




Gl 


0.43 (0.67) 


0.18 (0.56) 


0.12 (0.77) 


Rates+SK Spect 


G2 


0.43 (0.67) 


0.54 (1.0) 


0.83 (1.0) 




G3 


0.45 (0.72) 


- (0.57) 


- (0.93) 




Gl 


0.48 (0.78) 


0.28 (0.70) 


0.22 (0.83) 


Global fit 


G2 


0.48 (0.78) 


0.61 (1.0) 


0.83 (1.0) 




G3 


0.48 (0.83) 


0.20 (0.76) 


- (0.95) 



TABLE III. Maximum allowed value of COS23 cos^^ at 90 % and 99 % GL for the different 
solutions to the solar neutrino problem with the different statistical criteria. SMA is allowed at 90 
% for all the range of COS23 COS24. 
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FIG. 1. Allowed regions in Am^]^ and sin^(2t?i2) for the MSW four-neutrino oscillations 
from the measurements of the total event rates at Chlorine, Gallium, Kamiokande and Su- 
per-Kamiokande (825-day data sample). The different panels represent the allowed regions at 
99% (darker) and 90% CL (lighter) obtained as sections for fixed values of the mixing angles cf^c^^ 
of the three-dimensional volume defined by — Xmin~^-'^^ (90%), 11.36 (99%). The best-fit point 
in the three parameter space is plotted star. 
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FIG. 2. Allowed regions in Amli and sin^(2t?i2) for the vacuum four-neutrino oscillations 
from the measurements of the total event rates at Chlorine, Gallium, Kamiokande and Su- 
per-Kamiokande (825-day data sample). The different panels represent the allowed regions at 
99% (darker) and 90% CL (lighter) obtained as sections for fixed values of the mixing angles cf^c^^ 
of the three-dimensional volume defined by ~ Xmin"^-^^ (90%), 11.36 (99%) where Xmin 
the MSW region. 
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FIG. 3. Allowed regions in Am^]^ and sin^ i!}i2 for the MSW four-neutrino oscillations from the 
measurements of the event rates and the Super-Kamiokande zenith angular dependence data. The 
shadowed area represents the excluded region at 99% CL from the zenith angular data. 
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FIG. 4. Same as Fig]^ but for the measurements of the event rates and the Super-Kamiokande 
recoil electron energy spectrum. The shadowed area represents the excluded region at 99% CL 
from the energy spectrum. 
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FIG. 5. Same as Fig.^ but for the measurement of the event rates and the Super-Kamiokande 
recoil electron energy spectrum. The shadowed area represents the excluded region at 99% CL 
from the energy spectrum. 
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FIG. 6. Results of the global analysis for the allowed regions in /^.tti^i and sin^ 
four-neutrino oscillations. The light (dark) regions are allowed at 90% CL (99% CL) 
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FIG. 7. Results of the global analysis for the allowed regions in Am2i and sin^(2t?i2) for the 
vacuum four-neutrino oscillations. The light (dark) regions are allowed at 90% CL (99% CL). 
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FIG. 8. Ax as a function of the mixing pctrameter for the different solutions SMA (fuh 

Une), LMA (dashed), LOW(dotted) and vacuum (dot-dashed) from the analysis of the rates only 
right panels and from the analysis of the full data set (right panels). Each row represents the value 
for the different statistical criteria CI, C2 and C3 as defined in the text. The dotted horizontal 
lines correspond to the 90%, 95 %, 99% CL limits for each criteria. 
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